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Abstract 

Let  R^,...,Rm  be  rectangles  on  the  plane  with  sides  parallel  to  the 

coordinate  axes.  An  algorithm  is  described  to  find  the  contour  of 
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F  *  R.  LI  ...  U  R  in  0(mlogm  +  plog(2m  /p))  time,  where  p  is  the  number  of 
i  m 
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edges  in  the  contour.  This  is  0(m  )  (optimal)  in  the  general  case,  and 
O(mlogm)  (optimal)  when  F  is  without  holes  (then  p  <  8m-4) . 
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1.  Introduction 

In  this  paper  we  consider  the  following  geometric  problem-  Let 

R. , . . . ,R  be  rectangles  in  the  plane  with  sides  parallel  to  the  coordinate 
1  m 

axes.  The  union  F  *  R.  U  ...UR  may  consist  of  one  or  more  disjoint 

1  xn 

connected  components,  and  each  component  may  or  may  not  have  internal  holes 

(note  that  a  hole  may  contain  in  its  interior  some  connected  components  of  F) 

The  contour  (boundary)  of  F  consists  of  a  collection  of  disjoint  cycles 

composed  of  (alternating)  vertical  and  horizontal  edges.  By  convention,  any 

edge  is  directed  in  such  a  way  that  we  have  the  figure  on  the  left  while 

traversing  the  edge;  this  is  equivalent  to  saying  that  a  cycle  is  oriented 

clockwise  if  it  is  the  boundary  of  a  hole,  and  counterclockwise  if  it  is  an 

external  boundary  of  a  connected  component.  Given  R^-.-.R^,  our  task  is  to 

find  the  contour  of  F  *  R.  U  ...UR  • 

l  m 

There  are  several  applications  which  involve  iso-oriented  rectangles,  or 
without  loss  of  generality,  rectangles  whose  sides  are  parallel  to  the 
coordinate  axes.  Suffice  it  to  mention  its  relevance  to  Very  Large  Scale 
Integration  [5],  to  geography  and  related  areas  [3J,  to  computer  graphics 
(hidden- line  problems,  shadows,  etc.),  and  two-dimensional  data  organization 
[6].  There  is  also  considerable  theoretical  interest  in  problems  related  to 
the  one  discussed  in  this  paper,  such  as  finding  the  measure  of  the  union  of 
rectangles  [2],  testing  a  set  of  rectangles  for  disjointness  and  reporting 
all  the  Intersections  [4],  etc. 

2.  Informal  Description  of  the  Algorithm 

We  Informally  illustrate  the  method,  with  the  aid  of  an  example.  The 
algorithm  consists  of  two  phases.  In  the  first  phase  we  find  the  sec  V  of 
vertical  edges  of  the  contour  (edges  1  through  10  in  Fig .  1) ;  in  the  second 


2 


Fig.  1.  An  instance  of  the  problem. 


phase  we  link  these  vertical  edges  by  means  of  horizontal  edges  to  form  the 
oriented  cycles  of  the  contour. 

In  order  to  obtain  the  set  V  we  scan  the  abscissae  corresponding  to 
vertical  sides  of  rectangles  from  left  to  right.  At  a  generic  abscissa  c, 
the  section  J  of  the  vertical  line  x  ■  c  with  F  is  a  disjoint  union  of 
intervals.  This  section  remains  constant  between  two  consecutive  vertical 
sides  of  the  rectangles,  and  is  updated  in  the  scan  each  time  one  such 
side  is  reached.  If  s  is  a  vertical  side  of  some  rectangle  R  at  abscissa  c, 
the  portion  of  the  contour  of  F  contributed  by  s  is  given  by  s  H  J  where 
J  is,  as  usual,  the  union  of  intervals  representing  the  set  theoretical 
complement  of  J  on  the  line  x  •  c  (strictly  speaking,  J>  should  be 
understood  as  the  section  to  the  left  or  right  of  x  •  c,  depending  on 
whether  s  is  a  left  or  right  side  of  a  rectangle).  Storing  and  updating 
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of  J,  and  an  efficient  determination  of  s  H  J  require  the  use  of 
some  nontrivial  techniques  and  data  structures  and  is  deferred  to 
the  next  section.  We  assume  that  each  of  the  edges  in  V  is  directed 
upwards  or  downwards  depending  on  whether  it  originates  from  a  left  or 
right  side  of  a  rectangle. 

We  denote  by  (x;b,t)  a  vertical  edge  of  abscissa  x  having  b  and  t 
(b  <  t)  as  its  bottom  and  top  ordinates,  respectively;  similarly, 

(y;X,r)  represents  a  horizontal  edge  of  ordinate  y  having  l  and  r  (X  <  r) 
as  its  left  and  right  abscissae,  respectively. 

We  shall  now  describe  the  second  phase  of  the  algorithm.  We  first 
augment  the  set  V  by  adding  triples  (x;t,b)  for  all  (x;b,t)  in  V.  Notice 
that  every  triple  (x;y,,y2)  in  the  resulting  set  V*  is  in  a  one-to-one 
correspondence  with  a  vertex  (x,y^)  of  the  contour,  and  that  every  vertical 
edge  (x;b,t)  of  the  contour  is  represented  in  V*  exactly  twice  by  triples 
corresponding  to  vertices  (x,b)  and  (x,t)  of  the  contour.  We  now  sort 
the  triples  (x;y^,y2)  in  V*  in  lexicographic  ascending  order  by  (y^,x); 
letting  e^  ■  (x^b^.tj),  in  our  example  we  obtain  the  sequence: 

^X2’^2’C2^  ^x4*^4»^4^  (x9»^9»c9)  ^xio* t10’^310^ 

^x5 ’^5 ’ (X7*^7,ty)  ^xi’bi» tj_)  (x2»t2>l,2^  (x-^jt^jb^) 

(Xj^tpbp  (x3;b3,t3)  (Xgjbg.tg)  (Xgjbg.tg)  (Xgjtg.bg)  (xio*clobl(P 

(Xgjtg.bg)  (Xg.tg.bg)  . 

In  general,  let  a^,a2,...,ap  (p  even)  be  the  resulting  sequence.  It 
is  easy  to  see  that  the  horizontal  edges  of  the  contour  can  now  be  obtained 
as  the  segments  joining  the  vertices  repre  .ted  by  a2k_^  and  a2k*  ^or 
k  -  1,2,..., p/2.  More  exactly,  let  a2fc-1  -  (ijy.y^,  a2k  -  (r;y,y2) . 


T 


The  horizontal  edge  (y;i,r)  is  assigned  the  direction  from  left 
to  right  if  either  the  edge  corresponding  to  the  triple  (A; y,y^)  is 
directed  downwards  and  y  <  y^,  or  if  (A;y,y^)  is  directed  upwards  and 
y1  <  y;  otherwise  we  direct  (y;A,r)  from  right  to  left.  In  our  example, 
the  pair  (a^a^  *  ((x2;b2,  t2) ,  (x4;t>4,  t4> ) ,  with  b2  =  b4,  gives  rise 
to  the  horizontal  edge  (b2;x2>x4>  which  is  oriented  from  left  to  right 
because  the  edge  e2  -  corresponding  to  (x2;b2>t2)  -  is  directed  downwards 
and  b2  <  t2;  by  contrast  the  pair  (a17>a18)  “((XgJ tg.bg) ,  (x1Q;  t10>b1()) ) , 
with  tg  *  t^Q  gives  rise  to  the  edge  (tg;Xg,x^g),  which  is  directed  from 
right  to  left  because  eg  is  directed  downwards  but  tg  •<  bg.  It  is  clear 
that  by  a  single  scan  of  the  sequence  a^,...,a^  we  may  produce  all  the 
horizontal  edges  and  doubly  link  them  into  the  cycles  of  the  contour, 
in  such  a  way  that  the  ''forward"  links  determining  the  orientation  of  a 
cycle  are  consistent  with  the  direction  of  edges.  Then  we  may  find 
explicitly  the  cycles  of  the  contour  and  identify  each  of  them  by  a 
vertical  edge  with  the  minimal  value  of  abscissa  (this  can  be  done  in 
0(p)  time).  Notice  that  a  cycle  corresponds  to  the  boundary  of  a 
hole  if  this  edge  is  directed  upwards,  and  to  an  external  portion  of 
boundary  otherwise. 

3.  Efficient  Determination  of  s  Q  J 

In  order  to  efficiently  implement  the  first  phase  of  the  algorithm 
we  propose  to  normalize  the  ordinates  of  horizontal  edges  of  rectangles 
(i.e.,  to  replace  each  ordinate  by  its  rank  in  the  set  of  ordinates)  and 
to  make  use  of  a  segment  tree,  a  data  structure  introduced  by  Bentley  [2], 


which  we  now  recall. 
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Let  a,b  be  integers  with  a  <  b.  The  segment  tree  T(a,b)  is  built 
in  the  following  recursive  way.  It  consists  of  the  root  v  with  B[v]  *  a, 

E[v]  =  b,  and,  if  b-a  >  1,  of  a  left  subtree  T(a,  L(a+b) / 2 J )  and  right 
subtree  T(L(a+b)/2J  ,b) .  The  roots  of  these  subtrees  are  pointed  to  by 
LSONfv]  and  RSON[v],  respectively.  If  b-a  =  1  then  LSON[v]  =  RSON[v]  -  A. 

In  our  application  we  shall  use  the  segment  tree  T(l,h)  =  T,  where  h  i  2m 
is  the  number  of  distinct  ordinates. 

As  is  well-known  (see  also  the  procedure  INSERT  below),  an  interval 
s  =  [b,e]  (with  1  <  b  <  e  ^  h) ,  can  be  subdivided  into  O(logm)  segments,  each  of 
which  equals  [B[v] ,E(v] ] ,  for  some  node  v  of  T.  Thus,  for  each  v  we 
shall  have  an  integer  parameter  C[v],  which  denotes  the  current  multiplicity 
of  insertions  of  segment  [B[v],E[v]]  at  node  v.  We  shall  also  use  the 
parameter  STATUS [v],  which  may  assume  one  of  three  values:  full,  partial, 
empty.  Intuitively,  full  denotes  that  C(v]  >  0,  partial  denotes  that 
C[v)  *  0  but  C[u]  >  0  for  some  u  ^  v  in  the  subtree  rooted  at  v,  and 
empty  that  C[u]  *  0  for  every  u  in  this  subtree.  In  our  algorithm  we  scan 
the  vertical  sides  of  the  rectangles  from  left  to  right.  Each  left  side 
is  inserted  into  T  and  each  right  side  is  deleted  from  T,  so  that  the 
current  section  J  is  always  given  by  the  union  of  segments  [Btv],E[v]] 
over  all  full  nodes . 

Suppose  now  that,  for  some  node  v  in  T,  we  have  s  *  [b,e]  3  [B[v],E[v]]. 

The  contribution  of  v  to  s  D  J,  denoted  as  compl(v)  is  given  by 

f  0  if  STATUS  l v] -full 


compl(v)  -  [B(v],E[v]]  fl  J> 


compl(LSON[v])  U  compl(RS0N[v])  if  STATUS (v]-parti a 1 
lB(v],E[v]l  if  STATUS [v]-empty. 
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It  appears  therefore  that  s  fl  J  can  be  obtained  as  U  compl(v), 

v  6  U 

where  U  is  the  set  of  O(logm)  nodes  of  T  corresponding  to  the 
segmentation  of  s.  Notice  that  s  fl  J  may  consist  of  several  disjoint 
segments,  and  that  any  such  segment  is  obtained  as  a  collection  of 
contiguous  intervals  (corresponding  to  a  collection  of  0 (logm)  empty 
nodes  of  T) .  To  insure  that  each  contour  edge  be  produced  as  a  single 
item,  contiguous  intervals  have  to  be  merged.  To  implement  this, 
the  edges  of  s  fl  J  are  assembled  in  a  STACK  corresponding  to  a  bottom- 
to-top  scan.  At  the  top  of  STACK  there  is  always  the  upper  extreme  of 
the  last  inserted  edge  (or  initial  portion  of  this  edge);  if  the  lower 
extreme  of  the  next  segment  to  be  added  to  STACK  coincides  with  the 
top  of  STACK,  then  the  latter  is  deleted  from  STACK  prior  to  adding 
to  it  the  upper  extreme  of  the  next  segment,  thereby  realizing  the 
desired  merges . 

In  sunsnary,  an  interval  s  ■  [b,e]  is  inserted  into  T  by  the 
operation  INSERT (b,e, root (T) ) ,  where  INSERT(b,e,v)  is  the  following 
recursive  procedure: 

1  procedure  INSERT (b , e , v) 

(♦insert  segment  [b,e]  at  node  v  of  segment  tree*) 

2  begin 

3  if  (b  <  B(v])  and  (E[v]  ^  e)  then 

4  begin  COMPL(v)  (*see  below*) 

5  C[v]  :«  C[v]  +  1 

6  end 

7  else  begin  if  b  <  l(B[v]  +  E[v])/2J  then  INSERT (b,e,LSON[v]) 

8  L(B[v]  +  E[v])/2J  <  e  then  INSERT (b,e,RSON[v]) 

9  end 

(♦STATUS [v]  is  to  be  updated*) 

10  if  C[v]  >  0  then  STATUS [v]  :«  full 

11  else  if  LS0N(v]  ■  A  then  STATUS [v]  :*  empty  (*v  is  a  leaf*) 

12  else  if  STATUS [LS0N[v]]  -  STATUS [RS0N[v] 1  -  empty  then  STATUS [v]  *  empty 

13  else  STATUS [v]  :■  partial 

14  end 


7 


We  finally  give  Che  procedure  COMPL(v): 

1  procedure  COMPL(v) 

(*this  procedure  makes  use  of  a  STACK,  which  is  external  co  ic; 
it  pushes  on  STACK  a  sequence  of  segments  representing  the  set 
compl(v)*) 

2  begin 

3  if  STATUS [v]  =  empty  then  (*compl(v)  =  [B[v] ,E{v] ]*) 

4  begin  if  B[vl  =  top (STACK)  then  delete  top (STACK)  ^contiguous 

segments  are  merged*) 

5  else  STACK  *  B[v]  (*beginning  of  edge*) 

6  STACK  *  E[v]  (*current  termination  of  edge*) 

7  end 

8  else  if  STATUS [v]  =  partial  then 

9  begin  COMPL(LSON(v]) 

10  C0MPL(RS0N[v]) 

11  end 

12  end 

As  mentioned  before,  a  collection  of  vertical  contour  edges  is 
generated  either  in  correspondence  to  a  rectangle  left  side  (to  be 
inserted  into  T)  or  to  a  right  side  (to  be  deleted  from  T) .  The  procedure 
DELETE(e,b,v)  is  analogous  to  INSERT (e ,b ,v) .  The  only  difference  is  that 
in  lines  4  and  5  of  DELETE  we  have  C[v]  :=  C[v]-1  and  COMPL(v),  respectively 
(reversal  of  the  order  of  updating  C[v]  and  calling  COMPL(v)  reflects 
the  fact  that  for  a  right  side  the  J  in  s  0  J  is  understood  as  the  section 
to  the  right  of  the  current  abscissa) . 

Notice  that  in  order  to  guarantee  the  correctness  of  our  algorithm 
we  should  process,  for  any  value  of  x,  first  all  left  sides  (x;b,t)  in 
increasing  order  of  b,  and  then  all  right  sides  in  increasing  order  of  b. 

4.  Performance  Analysis 

We  shall  now  give  a  concise  description  of  the  algorithm.  Each 
is  represented  by  a  quadruple  Ci^r^jb^.t^),  where  b^  <  t^  and 

Ri  *  x 

Step  1:  Form  the  sets  (without  repetitions) 

A  »  [iL:  U  iS«}U  C  r  £ :  1  S  i  £  m},  B  -  [bt:  1515i]u  C^:  1515a], 
and  sort  each  of  them.  From  now  on  we  shall  identify  any  l ^  or  r^  with 
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its  rank  in  the  first  ordering,  and  any  or  t^  with  its  rank  in  the 
second  ordering.  (Comment :  In  this  way  we  transformed  our  problem  into 
an  equivalent  one  with  f ^ ,  r^  €  [1,...,|a|},  b^ , t^  €  [l, . . . , |B | }  for  1  ^  i  — 
Notice  that  |a(-,(b|  <  2m,  and  that  all  sorting  operations  are  now  doable 
in  0(m)  time  by  standard  bucket  sorting  (see  e.g.  [1]).  The  complexity  of 
this  step  is  clearly  O(mlogm)). 

Step  2:  Construct  tne  segment  tree  T.  (Comment:  This  requires  0(m)  steps.) 
Step  3:  Sort  the  set  (representing  all  vertical  sides  of  rectangles) 

S  =  [(Xi,left,bi,ti)  :  1  <  i  <  m)  U  [(r^ right, bi t^)  :  1  <  i  <  m] 
lexicographically  by  all  four  coordinates  (assuming  left  <  right)  and 
put  the  sorted  sequence  in  QUEUE.  (Comment:  This  can  be  done  in  0(m) 
time.) 

Step  4:  Scan  the  vertical  sides  of  rectangles  listed  in  QUEUE.  Any 
such  side  is  inserted  (by  calling  INSERT)  into  the  segment  tree  T  if 
it  is  a  left  side,  or  deleted  (by  calling  DELETE)  from  T  if  it  is  a 
right  side.  Each  of  these  op-  .ations  contributes  a  sequence  of  edges, 
placed  in  STACK,  which  has  to  be  emptied  (outputted)  each  time  we  pass  either 
to  the  next  value  of  abscissa,  or  from  left  to  right  sides  with  the 
same  abscissa  (this  is  necessary  in  order  to  avoid  incorrect  merging 
of  edges).  (Comment :  The  performance  analysis  of  this  step  will  be 
done  below.) 

Step  5:  Find  the  contour  cycles.  (Comment:  This  part  of  the  algorithm 
has  been  described  in  detail  in  Section  2  and  shown  to  run  in  time  0(p), 
where  p  is  the  number  of  contour  edges.) 

To  analyze  Step  4  we  shall  restrict  ourselves  to  calls  of  INSERT; 
calls  of  DEI£TE  can  be  handled  in  a  similar  way.  The  work  done  by  INSERT 
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is  best  analyzed  by  recalling  an  interesting  property  of  segment  trees. 

For  any  segment  [b,e],  let  i(b,e)  be  the  node  u  of  T  with  B[u]  =  b  and 

with  the  maximum  possible  E[u]  ^  e  (i.e.,  [B[u],E[u]]  is  the  left-most 

piece  in  the  segmentation  of  [b,e]);  similarly,  let  r(b,e)  be  the  node  u 

with  E[u]  =  e  and  the  minimum  possible  B[u]  £  b.  Denote  by  P  and  P  the 

Xj  it 

(sets  of  nodes  of  the)  paths  from  the  root  to  £(b,e)  and  r(b,e),  respectively. 

It  is  easily  seen  that  the  segmentation  of  [b,e]  is  given  by  i(b,e), 

r(b,e),  and  the  right  sons  of  P.\P  not  in  P,  as  well  as  the  left  sons 

h  r  Si 

of  p  \P  not  in  P  .  Since  T  is  a  nearly  balanced  tree,  the  length  of  each 
X*  IT 

of  its  paths  is  bounded  by  Tlogml  ;  thus  the  work  done  by  a  single  call  of 

INSERT(b,e,root(T))  -  disregarding  the  calls  of  COMPL  -  is  O(logm),  since 

it  corresponds  to  traversing  P^  U  P^,  and  spending  a  fixed  amount  of 

work  at  each  node  and  possibly  at  its  sibling. 

Now  we  shall  extend  the  analysis  to  account  for  the  work  done  by  COMPL. 

Let  (x;b,t)  be  the  rectangle  side  being  inserted  into  the  segment  tree. 

Let  [b,  t]  0  J  »  [c1,c2]  U  ...  U  C  c2fc—  1*  c2k  ^  *  C1  <  c2  <  “*  <  C2k'  let 

V2j-1  =  -4(c2j-l’C2j)’  v2j  *  r(c2j-l,c2j)’  1  ~  i  ~  k>  atld  let  P(vi)  be 
the  path  from  to  the  root  of  T.  It  is  easy  to  see  that  what  the 

recursive  procedure  INSERT  does,  together  with  COMPL,  is  exactly  the 

2k 

preorder  traversal  of  the  subtree  with  node  set  U  P(v  ).  A  fixed 

i»l 

amount  of  work  is  done  at  each  node  visited,  and  possibly  at  its 

2k 

sibling,  so  that  the  total  work  is  proportional  to  |  U  P(v, )J.  In 

i-1 

order  to  estimate  the  latter  quantity  we  shall  need  the  following 


lemma,  proved  in  Appendix  A: 

i 
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Lemma  1.  Let  X  be  a  binary  tree  with  N  leaves.  Assume  that  T  is 
balanced  so  that  any  path  from  the  root  to  a  leaf  is  of  length  L^°gNJ  or 
TlogNl  .  Let  L  be  a  subset  of  nodes  of  T,  and  for  every  u  €  L  denote  by  P(u) 
the  set  of  nodes  on  the  path  from  u  to  the  root.  Then 

j  U  P (u) j  <  | L| log  . 

Let  us  denote  by  n^  the  number  of  disjoint  pieces  in  0  J),  where  s. 

is  the  i-th  rectangle  side  processed  in  the  segment  tree  (1  <  i  <  2m) . 

By  our  lemma  applied  to  the  segment*  tree  (note  that  each  of  the  disjoint 

segments  [b,e]  in  H  J  contributes  at  most  two  members  of  set  L,  namely 

£(b,e)  and  r(b,e),  the  work  involved  in  processing  can  be  upperbounded 

by  Cn. log  for  some  constant  C  (this  expression  is  interpreted  as  0  for 
i  nt 

n^  -  0) .  Summing  it  over  all  i,  the  global  work  in  Step  4  is  upperbounded  by 

i  16m  ^  _  ,  64m^ 

C  2  n  log  — C  p  log  — —  , 

i-1  1  ni  P 

where  use  has  been  made  of  the  fact  that  2^  does  not  exceed  p/2,  the 
number  of  vertical  edges  of  the  contour  (strictly  speaking,  in  the 
presence  of  overlapping  vertical  rectangle  sides  we  should  write 
2n^  <  p/2  +  2m;  however,  this  does  not  affect  the  asymptotic 

2 

complexity  of  the  algorithm).  It  is  easy  to  see  that  p  <  m  +  4m,  so 
2  2 

that  C  p  log(64m  /p)  »  0(p  log(2m  /p)).  Putting  this  together  with  the 
previous  comments  on  Steps  1,  2,  3  and  5  of  the  algorithm,  we  obtain 
Theorem  2 .  The  algorithm  correctly  finds  the  p-edge  contour  of  a  union 


of  m  rectangles  in  time  0(mlogm  +  p  log(2m  /p)). 

2  2 
Notice  that  p  log (2m  /p)  is  increasing  in  p  in  the  range  [4,m  +  4m], 

for  m  >  4.  This  means  that  if  m  is  the  only  measure  of  the  size  of  the  input 


'  "V  '^5. 

*  r  '* 

'  T  .  ■ 
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then  the  complexity  of  our  algorithm  can  be  written  as 
2  2  2  2 

O(mlogm  +  (m  +  4m)log(2m  /(m  +  4m)))  =  0(m  ),  which  is  clearly  optimal 

2 

(the  algorithm  for  the  same  problem  in  [6]  is  of  complexity  0(m  logm)). 

Suppose  now  that  F  has  no  holes.  For  this  case  we  prove  in 
Appendix  B  the  following  lemma: 

Lemma  3.  Let  F  =  R^  U  . . .  1)  be  without  holes,  and  suppose  it  is 

composed  of  k  disjoint  connected  components.  Denote  by  p  the  number  of 

edges  in  the  contour  of  F.  Then  -p  ^  8m-4k. 

Therefore,  in  the  absence  of  holes,  the  algorithm  runs  in  time 
2 

O(mlogm  +  8mlog(2m  /p))  *  O(mlogm).  We  shall  now  prove  that  this  is 
optimal,  using  as  a  computation  model  the  RAM  extended  with  real- number 
arithmetic. 

Theorem  4.  The  complexity  of  finding  a  contour  of  F  *  R.  U  ...UR  , 

. .  1  m 

where  F  is  without  holes,  is  Q(mlogm). 

Proof.  We  shall  show  that  sorting  of  m  numbers  x, ,...,x  reduces  in 

0 (m)  time  to  our  problem.  Indeed,  given  x^,  let  R^  be  the  cartesian 

product  of  the  x-interval  [0,x^]  and  of  the  y-interval  [0,M-x^],  where 

M  ■  maxj  <  ^  <  mx^  +  1  (see  Fig.  2;  without  loss  of  generality  we  assume 

x,,...,x  >0).  It  is  clear  that  F  *  R.  U  ...  U  R  is  without  holes  and 
in  i  m 

that  from  the  contour  of  F  we  can  obtain  the  sorted  sequence  of  the  x^  in 
0 (m)  time .  O 


aw 
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IC  is  easy  to  see  that  TQ(L)  is  a  tree  with  |l|  leaves.  It  is  also 
clear  that  |T(L)|  is  maximum  (for  a  fixed  |l|)  if  TQ(L)  is  a  balanced 
top-most  part  of  T,  with  |l|  leaves  and  2 j l| -1  nodes.  In  such  an 
extremal  case 

jT(L)  |  <  2 1 L| -1  +  ]  L|  ( T log  Nl  -  Llog  |  l|j  ) 

<  |L| (2  +  log  N  +  1  -  log  J  L |  +  1)  =  j  l| log  .  □ 


rectangles  are  labelled  so  that  i.  ^  L  ,  for  1  <  i  <  m.  It  is  easy  to 

i  mrl 

see  that  F  *  E.  U  ...IIR  is  without  holes.  Indeed  R  cannot  fill  a  hole 
1  m  tnri 

in  F,  because  lies  entirely  to  the  right  of  any  hole  in  F.  Suppose  F 

consists  of  k  connected  parts,  and  R  ^  "glues  together"  q  parts  of  F. 

Notice  that  it  necessarily  means  that  the  intersection  of  the  left  side 


of  R^^  with  F  consists  of  q  disjoint  segments.  It  is  easy  to  see  that 
the  increase  of  the  number  of  edges  in  the  contour  caused  by  adding 


Rnrt-1 


may  come  from  the  following  sources  (see  Fig.  4). 
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(a)  segments  between  the  q  disjoint  segments 

on  the  left  side  of  R  : 

m+1 

(b)  similar  segments  on  the  right  side  of  R^^: 

(c)  bottom  and  top  side  of  R^^: 

(d)  doubling  the  edges  of  the  contour  of  F  which 
intersect  both  the  left  and  right  side  of  R^^: 

TOTAL  at  most  4q+4 

Taking  into  account  that  the  number  of  connected  components  of  F* 
is  equal  to  k*  »  k-q+1,  and  using  the  inductive  assumption,  we  obtain  the 
following  estimate  on  the  number  of  edges  in  the  contour  of  F: 
p*  <  (8m  -  4k)  +  (4q  +  4)  »  8(m  +  1)  -  4(k  -  q  +  1)  *  8(m  +  1)  -  4k*. 


at  most  q+1 
at  most  q+1 
ac  most  2 

at  most  2q. 


Notice  that  there  exist  hole-free  unions  of  m  rectangles  with 
p  -  8m  -  4k,  for  any  m  and  any  1  5  k  <  n.  Every  connected  part  of  such 
a  union  is  of  the  form  shown  in  Fig.  5. 


Fig.  5.  A  hole-free  connected  union  of  rectangles  with  p  *  8m  -  4. 
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